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Abstract. In this letter we present first results from spec- 
troscopic observations of Diffuse Ionized Gas (DIG) in 
the halo of NGC 5775 obtained with FORS1 1 attached to 
UTl/Antu of the Very Large Telescope (VLT). At our slit 
position perpendicular to the disk (41" SE of the nucleus) 
the emission of [Nn]A6583, [O m]A5007, and [O n]A3727 is 
detected out to 9kpc into the halo, allowing possible ion- 
ization mechanisms of the DIG to be examined. Photoion- 
ization models which assume a dilute radiation field are 
able to fit the data for the disk, but they cannot account 
for the line ratios measured in the halo (e.g., [Ol]/Ha or 
[On] /Ha). In particular they fail to predict the observed 
increase of [O in] A5007/Ha with increasing \z\. The most 
striking result concerns the kinematics of the halo gas. Ve- 
locities at high galactic latitudes drop from the midplane 
value to reach the systemic velocity at z m 9 kpc. An anal- 
ysis of VLA archive data of the polarized radio-continuum 
emission at 4.86 GHz and 1.49 GHz reveals that magnetic 
fields in the halo have a strong component perpendicu- 
lar to the disk and are aligned with the Ha and radio- 
continuum spurs in the halo. This can result either from 
a strong wind action or, more likely, from the generation 
of dipolar magnetic fields. We briefly discuss the interre- 
lation of the magnetic field structure and gas dynamics, 
in particular the role of magnetic fields in gas outflows, as 
well as the possible heating of DIG by magnetic reconnec- 
tion. 

Key words: Galaxies: halos - Galaxies: individual 
(NGC 5775) - Galaxies: ISM - Galaxies: kinematics and 
dynamics - Galaxies: magnetic fields - Galaxies: structure 



1. Introduction 

There is increasing evidence that halos of normal spiral 
galaxies contain a complex multiphase interstellar medium 
(ISM) (for a review see, e.g., Dahlem 1997 ). Besides a 
hot component as observed by X-ray emission (>10 6 K) or 
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UV absorption (~0.5xl0 5 K), Diffuse Ionized Gas (DIG) 
constitute s ano ther important phase of typically 10 4 K 
(Dettmar 1992 ). This component is easily observed by its 
H II region like line emission and is found to be wides pread 
in galaxi es wi th actively star for ming disks (Rand 
Dettmar |1999| , Rossa k Dettmar fcOOOp . 



1996 



In the last couple of years two important problems 
have been discussed with regard to the physics of this 
interstellar gas phase. The observed line ratios are not 
explained in a straight forward way by pure photoion- 
ization models and require additional heating mcchan- 
sims of the halo gas (e.g., Golla et al. [19961 Rand |1997, 



Tiillmann & Dettmar 2000D while first kinematic observa- 



tions hint at a slower rotation of halo gas if compared to 
the underlying disk (Rand [1998 ). Both kinds of observa- 
tions could help to constrain models of the origin of halo 
gas which in turn are important for the understanding of 
galaxy evolution since large scale mass flows are of influ- 
ence for star formation and chemical evolution. 

Only for a few objects, such as NGC 891 (Dahlem et al. 



199J), NGC 4631 (Golla & Hummel |199J) , and NGC 5775, 
the various components of the halo ISM are observed with 
sufficient sensitivity and resolution to study possible phys- 
ical processes at their interfaces such as turbulent mixing 
or thermal conduction. For such processes yet another ISM 
component can be of importance, namely the non-thermal 
component consisting of magnetic fields and cosmic rays 
as observed by its radio-continuum radiation. Any large 
scale gas motion and in particular its z-component will 
also be of importance for dynamo theories which try to 
explain the large scale distribution of magnetic fields in 
galaxies. However, various approaches to describe galactic 
dynamos are still very controversial and more detailed in- 
formation on galactic gas flows in halos as well as magnetic 
field structures could help constraining such models. 

For all these questions, the case of the Sc edge-on 
galaxy NGC 5775 has received some attention recently, 



since marked radio-continuum (Duric et al. 1998) and Hi 



(Irwin 1994) halo components are accompanied by a very 
extended DIG halo (Dettmar [1992[ Collins et al. |2000p 
for which a strong decline of the rotational velocity with 
height above the plane (z) has been found (Rand 2000). 
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In this letter we report on first spectroscopic results 
for the DIG halo of NGC 5775 with the VLT, confirming 
the strong z-dcpcndcncc of the halo rotation presented by 
Rand ( |2000D , along with a first detailed determination of 
the magnetic field structure in the halo from VLA radio- 
continuum observations. 



2. Observations 

Spectroscopic observations were carried out with FORSl[] 
at VLT/UT1 in longslit mode on 1999 June 14-15 at a me- 
dian seeing of 0'.'7. A Tek 2048 x 2046 CCD with a pixel size 
of 24 x 24/xm and a spatial resolution of C/2 pix - 1 was used 
for the measurements. Grism GRIS_600B+12 was chosen 
for the blue and GRIS_600R+14 for the red spectral re- 
gion, yielding a wavelength coverage of 3500-5900A and 
5250-7400A, respectively. The slit length of FORS1 is 6f8 
and the selected slit width of 1'.'31 gives an average resolu- 
tion of 5.9A for the red and 6.2A for the blue grism, respec- 
tively. Our slit covers large filamentary structures of Ha 
emission in the NE and SW halo (Fig. Q). For each grism 
the total integration time amounted to 3 hours, reaching a 
noise level of 1.26 x 10 -19 ergs cm -2 s _1 A -1 pix -1 . Data 
reduction was performed in the usual manner using stan- 
dard IRAF software while the analysis was done within 
the MIDAS environment. 

Maps of total power and the linearly polarized 
radio-continuum emission of NGC 5775 at 1.49 GHz and 
4.86 GHz have been obtained from archive data of the Very 
Large Array (VLA) of the National Radio Astronomy Ob- 
servatory (NRACQ). All data were recorded during 1989 - 
1994. The total observation time used on source was about 
10 hours at 4.86GHz (D-array), 14hours at 1.49GHz with 
C-array, and 2.5 hours with D-array. Data reduction has 
been performed with the AIPS software package following 
standard procedures. The r.m.s. noise obtained in the fi- 
nal maps is 16/iJy/b.a. and 100/iJy/b.a. in total power, 
and 7/xJy/b.a. and 25/iJy/b.a. in polarized intensity at 
4.86 GHz and 1.49 GHz, respectively. 



3. Results and discussions 

3.1. Ionization 

Several strong emission lines such as [Nn]A6583 or 
[On]A3727 can be traced out to 66" or 9kpc (assuming 
a distance of D = 26.7 Mpc) above the disk in both di- 
rections, much further than any previous imaging result. 
Similar findings have been reported by Rand 2000 for the 
western part of NGC 5775. All measured line ratios are 
presented in Fig. ||. Some ratios, such as [S 11] and [Oi], 
are lower in the SW section, most likely due to a more 
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Fig. 1. Continuum subtracted Ha image of NGC 5775 ob- 
tained with EMMI at the NTT. Ha-contours shown for 
the disk reveal that si covers a gap of Ha emission in the 
disk, cutting also through the NE and SW filament. 

patchy structure of the filament. With an inclination of 
i = 84°, NGC 5775 is close to edge-on and inclination ef- 
fects are considered to be negligible. In order to clarify 
the issue of possible ionization mechanisms of the DIG we 
compare our observational results with predictions of pure 
photoionization models from Mathis (1986, Ma86 here- 
after) and Domgorgenfe Matins ( |l994| , DM94 hereafter). 
If we adopt photoionization by OB stars as the only ioniza- 
tion source and that radiation gets more and more dilute 
as it extends towards the halo, these models should repro- 
duce observational data, such as the rising [O ill] /Ha, non- 
linear changes in [S n]/[N 11], or extreme ratios of [O i]/Ha 
and He i/Ha. 

[Nil] /Ha is fitted only for the disk by the model of 
Ma86 assuming a very soft radiation field (log q — —6) and 
07 stars with effective temperatures of 35000 K as ioniz- 
ing sources, whereas the DM94 model fails completely. By 
no means DM94 can reproduce halo ratios for [Sii]/Ha 
larger than 0.66, while Ma86 gives correct values ranging 
from 0.26 for the disk to 0.88 for the halo. Since these 
models predict an increase of [S 11] /[N 11] towards the halo, 
most likely a consequence of slightly different ionization 
potentials of S + and N + , they are unable to explain the 
non-linear trend visible in Fig. || Another important ra- 
tio is Hei/Ha which could be measured only close to the 
disk-halo interface of NGC 5775 (-1.1 kpc < z < 1.65 kpc), 
making this detection nevertheless unique concerning its 
z extent. However, both models underestimate this ratio 
by a factor of two. 

The most striking feature in Fig. |^ is the rise of 
[O ill] /Ha with increasing z which is not explained by pure 



R. Tiillmann et al.: The thermal and non-thermal gaseous halo of NGC5775 



3 



0.8 

o.e 

0.4 
O.S 
0.8 
0.6 
0.4 
OS 
O 
1 

0.8 
0.6 
0.4 

0.15 
O 1 

0.05 



i - 



I I I I I I I I I I I I I I I I I 



I I I I I I I I I I I I I I I I I I I 



T gj »«„-"» . £ 

1 *S %r » 



I I I I I I I I I I I I I I I I I I I 



- lO —8 —6 —4 — 2 O 2 4 6 8 lO 
SW z (kpc) NE 



0-8 

0-7 

0.6 

O.S 
1 

0.8 
0-6 
0-4 
0-2 
O 

o.e 

o.e 

0.4 
0.2 

o 

O.S 
0.4 
0.3 
0.2 
0. 1 







J 1 1 1 1 1 1 1 1 1 1 1 1 

- O.H5 .... . ™ ' ot^'- 

H >K 1 ' 1 


1 1 1 L 

.¥ '"$.§ / - 

1 h "1 1 ; 


7 i ?*i It 
-o.so %^ /••- 

; 1 1 1 1 1 1 1 1 ryy 1 . 1 


.... 1 1 1 1 1 ; 


- s 





-IO —8 —6 —4 —2 O 2 4 6 8 1C 
SW z (kpc) NE 



Fig. 2. Line ratios (corrected for interstellar extinction) plotted as a function of \z\. Theoretical line ratios (dotted 
lines) based on temperature profiles (measured from nitrogen emission, upper right) are fitted to the data by variing the 
corresponding ionization fraction. "x"-symbols are from temperature determinations using upper limits of [Nn]A5755 
emission and are given as a consistency check. 



photoionization codes. Generally this ratio declines from 
values of about 2 in the disk to 0.8 in the halo. 

Both models also fail to predict [On] /Ha. The lowest 
predicted values vary around 1, whereas the observed data 
never exceed this limit. [Oi]/Ha reveals with 0.33 for the 
SW and 0.48 for the NE halo by far the highest values 
ever measured for the DIG. Ma86 and DM94 are able to 
fit ratios only up to 0.1. 

The extreme line ratios of Hei/Ha (>0.05) and 
[Oi]/Ha (>0.3), the shapes of [Om]/Ha and [Sii]/[Nii] 
together with an unprecedented DIG extent of 9 kpc (see 
also Rand 200C| ) make this galaxy the most outstanding 
and challenging testbench for future ionization models. 

Since common photoionization models are unable to 
reproduce the data for halo DIG in external galaxies as 
well as the Reynolds-layer of the Milky Way, recent stud- 
ies (e.g., Rand |l998j, Reynolds et al. |l999j, Tiillmann & 



Dettmar 200C) have pointed out the need to involve ad- 
ditional ionization and/or heating mechanisms such as 
shocks, photoelectric heating by dust, or magnetic recon- 
nection. These mechanisms should increase the electron 
temperature without affecting the ionization stage of an 
atom and dominate in the halo over collisional ionization. 
Following Haffner et al. fll999| ) and Collins & Rand ( |2000 , 
CR hereafter) we used the [Nn]A6583 emission line to de- 
termine a temperature profile (shown in the upper right 
diagram of Fig. ||) which varies between 6500 and 9200 K. 
This profile can be used to derive theoretical line ratios 
(dotted lines in Fig. |^) for oxygen and sulfur. In order 
to fit the observed ratios, the ionization fraction as the 
only free parameter is changed accordingly. For consis- 



tency we also derived electron temperatures for the east- 
ern part of NGC 5775 which are based on upper limits 
of the [Nn]A5755 emission line (Osterbrock, |l989| ). They 
have been estimated to be 7350 K in the disk (at 0.5 kpc), 
7750 K for the NE filament, and 7900 K for the SW fila- 
ment (both measured at 5.0 kpc), respectively. These val- 
ues are in good agreement with our temperature determi- 
nations mentioned above and also with profiles obtained 
by CR for the central and western part of this galaxy. 

Furthermore, plots on the right of Fig. |^ reveal promi- 
nent gradients which decrease as the ionization stage of 
oxygen increases. A likely physical reason could be a ris- 
ing te mpera ture towards the halo, according to Reynolds 
et al. fll99Sp . This can simply be shown if we plot theoret- 
ical oxygen line ratios per ionization fraction vs. electron 
temperature. As a result [Oi]/Ha increases much faster 
with T than [Oil] /Ha or [Oin]/Ha. Hence, different tem- 
perature dependencies of oxygen ionization stages are able 
to reflect the observed gradients. 

3.2. Kinematics 



The most surprising result revealed by our deep spectra 
concerns the kinematics of DIG at large z. The bright 
emission lines of hydrogen, nitrogen, and oxygen have been 
used to derive heliocentric velocities for the DIG as a func- 
tion of z (vertical distance above the plane). As can be 
seen in Fig. || the gas velocity is with 1822 ± 8kms _1 
highest at z = 0kpc and —1.6 kpc, decreasing slowly to- 
wards larger z distances, finally reaching values cl ose or 
equal to the systemic velocity of 1681km s _1 (Irwin 1994 ) 
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Fig. 3. Heliocentric velocities as a function of z for the 
minor slit si. Dashed lines indicate the systemic velocity 



as measured by Irwin (1994). No qualitative differ 



ences are visible between the kinematical behavior of Ha, 
[Om]A5007, and [Nn]A6583. 



at z = 9 kpc. The measured trend is explainable assum- 
ing the DIG not to rotate at high galactic z. A similar 
result was recently reported by Rand ( 2000| ) for two dif- 
ferent slit positions. This allows us to examine whether 
the halo DIG rotates regularly since our slit position and 
slit 1 of Rand (200C) both have a distance of 45" from 
the center corresponding to 6.1 kpc but on opposite sides 
of the disk. The maximum rotational velocities at these 



positions (with respect to v sys at z — 2 kpc) deviate only 
by 4kms _1 and the dependence of v ro t(z) is to first order 
in agreement with a projected regular velocity field. How- 
ever, different gradients imply that the kinematics of the 
DIG halo is more complex than a symmetric and static 
model would allow for. This requires more modelling in- 
cluding the shape of the dark matter potential and pos- 
sible influences of outflow kinematics. The resolution of 
our spectra unfortunately is unsufficient to check for the 
velocity dispersion of the lines. 

Since our spectra have very high sensitivity at very 
good spatial resolution, Fig. || also indicates some small- 
scale structure in the velocity field with localized minima 
at z — —8.5, —4.0, 1.5, 5.2, and 9.5 kpc. Close to the disk, 
in particular in the NE, these changes can be explained by 
dust obscuration, "shifting" the peak velocity at the cen- 
ter to z = — 1.6 kpc. An R-band image of NGC 5775 gives 
evidence for large amounts of dust, mainly in the northern 
part, extending from the disk plane 1.8 kpc into the halo. 
The kinematical structures at larger z are not directly as- 
sociated with observed dust features and hint at either an 
inhomogeneous DIG distribution along the line of sight or 
effects from flows. 



3.3. Magnetic field structure in the halo of NGC 5775 

In Fig. H we present maps of total power, polarized in- 
tensity, and apparent polarization B-vectors at 4.86 GHz 
overlaid upon an enhanced Ha image. NGC 5775 shows an 
extended radio envelope which can be traced beyond 2 kpc 
from the disk plane. While close to the disk plane the po- 
larization B-vectors are generally disk parallel, at heights 
above 1 kpc they form a X-like pattern, corresponding also 
to similarly shaped extensions of polarized intensity. 

Before discussing the magnetic field structure we make 
a comparison of the polarization information at both fre- 
quencies to estimate the possible influence of Faraday ro- 
tation and depolarization. For heights less than 1 kpc our 
estimates for a regular azimuthal magnetic field oriented 
parallel to the disk, assuming pressure balance with cos- 
mic rays and typical thermal electron densities of 0.03 
cm" 1 , imply that NGC 5775 is Faraday thick at 4.86 GHz 
close to the disk plane. Therefore B-vector orientations 
are not conclusive. However, a reasonable agreement be- 
tween the orientations of polarization angles at 4.86 GHz 
and 1.49 GHz, a relatively high polarization degree (reach- 
ing 30%) at 4.86 GHz, and the depolarization at 1.49 GHz 
by a factor of about 0.3-0.4 support the assumption that 
above 1 kpc the halo of NGC 5775 is Faraday-thin down 
to 1.49 GHz. Thus, if the B-vector orientations at two 
frequencies differ there locally by some tens of degrees, 
the rotation angle at 1.49 GHz is smaller than 90°, hence 
the Faraday rotation bias at 4.86 GHz does not exceed 
10°. However, in the SW part of the halo (R. A. 2000 — 
14 h 53 m 54 s Dec 20 oo = +03°33'15"), similar orientations 
of B-vectors at both frequencies are apparently in conflict 
with the depolarization at 1.49 GHz by a factor of 0.1, 
which raises the suspicion that it might be Faraday-thick. 
On the other hand numerical simulations show that this 
depolarization can also be due to pure Faraday dispersion 
in random fields, without substantial Faraday rotation. 

From the above comparison we conclude that the ob- 
served B-vectors at 4.86 GHz in the halo at heights \z\> 
1 kpc are not significantly affected by Faraday rotation 
and that therefore the magnetic field in the halo has a 
significant vertical component. Such a vertical magnetic 
field component strong enough to be seen in emission can 
be generated in two ways. A spherically blowing galac- 
tic wind could result in an enhanced vertical component 
of the magnetic field. While the poloidal magnetic field 
preserves its quadrupole symmetry, resulting in B-vectors 
parallel to the plane below z ~ 1 kpc, characteristic X- 
shaped structures develop at larger heights (Brandenburg 
et al. 



1992) 



An alternative way to obtain the vertical B-vectors is 
the generation of the dipolar (called A-type) poloidal mag- 
netic field produced by the dynamo process. The principal 
condition is a rigid rotation extending over a substantial 
range of galactocentric radius and a large vertical scale- 
height of the ionized gas, both being true for NGC 5775 
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Fig. 4. VLA radio-continuum maps with total power (left), polarized intensity (right) and resulting B-vectors overlaid 
on a Ha image of NGC5775. Contours are at 3, 8, 21, 55, 144, 377, and 610xl6^Jy/beam for TP and 3, 5, 8, 13, 21, 
and 34x7^Jy/beam for PI, for the B-vectors a length of 1" corresponds to lO^Jy. 



(Lehnert & Heckman 1996| ) . In this case the B-component 
parallel to the disk would vanish. Unfortunatley, this way 
of distigushing both mechanisms cannot be used here be- 
cause of the strong Faraday effects close to the disk. How- 
ever, we note that the detailed analysis of Faraday rotation 
angles reveals field reversals across the plane (excluding 
z < 1 kpc) in favour of a dipolar poloidal field. 

Whether a vertical decrease of rotational speed, as seen 
in Fig. ||, makes this field mode growing faster remains yet 
to be investigated but cannot be excluded. We note that 
the existing dynamo models which account for rotation 
speed falling with z yield very different results. While lit- 
tle effect of vertical rotation decrease upon the dynamo 
efficiency was obtained for a classical dynamo concept 
(Brandenburg et al. 1993) a special model of a supernova- 



NGC5775 TP 1.49GHz filtered 



vectors 4.86GHz 



driven dynamo by Ferriere & Schmitt ( 2000 ) fails in this 



case to generate any stable, growing magnetic field mode. 
The dynamo models in case of a vertical rotation drop-off 
can be validated if only more information on the magnetic 
field structure in halos along with kinematical information 
would be available. 

However, as noted by Collins et al. ( |2000 ) we need 
to consider the influence of the vertical magnetic field 
structure on the cosmic ray propagation and the gas 
flows. To clearly demonstrate the correlation of the radio- 
continuum spurs with the magnetic field structure we 
have applied a median filter to the total power map at 
1.49 GHz, which revealed structural components inclined 
to the galaxy plane in the same X-shaped manner as the 
polarization B-vectors (Fig. ||). 

The polarization features are anchored in the disk 
plane with the SE one being attached to a local radio 
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Fig. 5. The strongly filtered total power image of 
NGC5775 at 1.49 GHz with a median filter applied. It 
shows in detail the loci of brightness maxima in the dis- 
tribution of the total power brightness. 



bright region in the disk, most likely a star forming re- 
gion. The total power spur at R.A. 20 oo = 14 h 53 m 53!l, 
Dec2ooo = +03°32'55" curves towards south in its outer 
part; the same behaviour is seen for the polarization vec- 
tors at both frequencies. We suspect that they are associ- 
ated with cosmic ray electron streaming from intensively 
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star-forming regions along the inclined magnetic lines. The 
inclined spurs extend down to the disk plane which allows 
us to speculate that magnetic lines may do the same, as 
expected for the A-mode. We also note that the cosmic 
ray propagation along the magnetic lines of the regular 
field may yield an X-shaped distribution of the polarized 
intensity (Fig. ^). Similarly, an easier streaming of the 
ionized gas along magnetic fields of the discussed geome- 
try can give rise to the occurrence of Ha-emitting spurs 
coincident with the radio ones and to an increased verti- 
cal scale of ionized gas in regions of magnetic structures 
highly inclined to the disk plane. 

Finally we want to address the question whether the 
energy stored in the magnetic field could in principal help 
solving the heating problem discussed in sect. 3.1. Mag- 
netic reconnection as a he at sou rce for DIG was sug geste d 
earlier, e.g., by Birk et al. ( 1998| ) or Reynolds et al. ( 1999| ), 
the latter based on an analogy from solar physics pro- 
posed by Raymond (1992). In the following we use the 



heating rate for magnetic reconnection G mag as derived 



by Lesch & Bender ( |1990| ): G mag = B 2 v a /8ttL, with v A 
the Alfen velocity and L the dissipation length. For char- 
acteristic regions in the spurs the minimum total magnetic 
field strength implied by the pressure balance with cosmic 
rays is > 5 /iG. With this value we obtain the required 



heating rates as determined by Reynolds at al. (1999) for 
DIG densities of typically n c = 1CP 3 cm~ 3 if L is in the or- 
der of parsec, not unreasonable for structures in the ISM. 
Since this heating rate depends on the Alfen velocity and 

— 1/2 

thus on n e the increasing importance of the additional 
heating source with decreasing particle density would be a 
natural consequence. Reconnection could also be respon- 
sible for the growth of ordered magnetic field structures 
in the spurs. It has been suggested that this process could 
help to reorganize the field structure in a galactic fountain 



flow (e.g., Kahn 1991) 



4. Summary and conclusions 

The most important results can be summarized as follows: 

— emission lines of DIG can be traced out to 9 kpc in the 
halo of NGC 5775 

— the rotational velocity of the DIG drops from the rota- 
tional value at the midplane and approaches the sys- 
temic velocity at distances z > 6.0 kpc 

— pure photoionization models are unable to reproduce 
the extreme line ratios of [Oi]/Ha, Hei/Ha, and 
[O in] /Ha as measured for the halo 

— a high degree of polarization and a strong vertical com- 
ponent of the magnetic field are observed in the halo. 
We found some, yet weak evidence for the dipolar mag- 
netic field 

— a correlation exists between extraplanar DIG and mag- 
netic fields, evidenced by the alignment of magnetic B- 
vectors with prominent halo features in the DIG and 
radio-continuum distribution 



— we suggest that reconnection processes play a substan- 
tial role in both heating the DIG structures in the halo 
and building up regular fields away from the disk plane 
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